Product Wave Function 
Renormalization Group 



T. Nishino 1 ) and K. Okunishi 2 ) 

(1) Physics Department, Graduate School of Science, 

Tohoku University, Sendai 980-77 

(2) Physics Department, Graduate School of Science, 

Osaka University, Toyonaka, Osaka 560 



Synopsis 

The authors propose a fast numerical renormalization group method 
- the product wave function renormalization group (PWFRG) method - 
for ID quantum lattice models and 2D classical ones. A variational wave 
function, which is expressed by a matrix product, is improved through a 
self-consistent calculation. The new method has the same fixed point as the 
density matrix renormalization group method. 
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The real-space renormalization group 1 ) (RSRG) is a basic concept in 
statistical physics. One of the recent progress in the RSRG is the density ma- 
trix renormalization group (DMRG) method established by White. 2 ' 3 ) The 
DMRG method has been applied to various one-dimensional (ID) quantum 
lattice models, such as the spin chains 3 ) and electron systems, 4 ' 5 ) because the 
method gives precise results for large scale systems. The method is also ap- 
plied to 2D classical models. 6 ) Recently, fixed point structure of the DMRG 
method has been analyzed by Ostlund and Rommer. 7 ) They show that the 
ground state wave function obtained by the DMRG method is a product of 
matrices. 

Quite recently, the authors have shown that the DMRG method is 
essentially the same as Baxter's variational method on the corner trans- 
fer matrix. 8 ) On the basis of this fact, we have formulated a fast numeri- 
cal method — the corner transfer matrix renormalization group (CTMRG) 
method 9 ) — for 2D classical systems. The CMTRG method is applicable to 
ID quantum lattice models with the help of the Trotter formula. 10 ) 

In this paper, we propose a new renormalization group method, which 
is an extension of the CTMRG method. The new method is applicable 
to ID quantum systems with no reference of the Trotter formula. 10 ) Since 
the method renormalizes a wave function, which is expressed by a ma- 
trix product, 11 ) we call it 'product wave function renormalization group 
(PWFRG) method.' Since the PWFRG method and the DMRG method 
have many aspects in common, we review the DMRG method at first. We 
then show the numerical algorithm of the PWFRG method for 2D classi- 
cal models. We accelerate the numerical calculation with the help of the 
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modified Lanczos method. 12 ' 13 ) Finally we discuss the way how to apply the 
PWFRG method to ID quantum systems. 

We choose the 'interaction round a face (IRF) model' as an example of 
the 2D classical model. 8 -* The IRF model is defined by a Boltzmann weight 
W(a'b'\ab) on each face, that is surrounded by four n-state spins a, b, a', and 
b'. The row-to-row transfer matrix is 

T(a'b'c . . . y'z'\abc ...yz)= W(a'b'\ab) W(b'c'\bc) . . . W(y'z'\yz), (1) 

where the variables {a'b'c' . . .y'z'} and {abc. . .yz} denote the n-state spins 
in subsequent rows. We assume that W(ab\cd) is symmetric — W(ab\cd) = 
W(ba\dc) = W(ca\db) = W(dc\ba) in order to simplify the discussion. 
Generalizations for asymmetric cases are straightforward. 

The DMRG method maps the transfer matrix T in Eq.l into the effective 
one 6 ) 

fgi'j'r/ltijTi) = Pii'emWii'j'lijWv'tiv), (2) 

where P represents the effective transfer matrix for the left/right-half lattice. 
The greek indices £, rj and rj' denote m-state block-spin variables, that are 
shown by squares in Fig.l. The eigenvector V that corresponds to the largest 
eigenvalue of T is well approximated by a matrix product 2 ' 3 ' 7 -' (Fig.l) 

Vmjv) = 52R(i£\a)u; a R T (a\jTi), (3) 

a 

where R is a nm by m orthogonal matrix, R T is its transpose, and the 
relation R T (a\i^) R(i^\(3) = 5@ is satisfied. We have expressed the 
(nm) 2 -dimensional vector V, which we call 'product wave function,' as a 
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nm-dimensional square matrix. The r.h.s. of Eq.3 is actually the singular- 
value decomposition of the l.h.s of Eq.3, and therefore u a is the singular 
value (or the eigenvalue) of the 'matrix' V{i£\jrj). 

The PWFRG method gives the fixed point values of P and R — those 
in the thermodynamic limit — through successive improvements. We start 
from the 4-site system. 2 ' 3 ) The initial conditions for the system with open 
boundary conditions are P = W (that means T = WWW), R{ia\(3) = 5^5^, 
uja = 1/y/n, and m = n. We have to set appropriate initial values for P 
and u a when fixed boundary conditions are imposed. Starting from the 
initial status, we improve P and R by way of the following self-consistent 
calculations. 

As a first step, we create a trial product wave function by using R and 
<jj a : V{i£\jrj) = ^2 a R(i$ t \a)u> a R T (a\jr]). The vector V is not usually the 
eigenvector of T = PWP in Eq.2. We then multiply T on V in order to 
obtain the improved Ritz vector 

V'tfZW) = ^P(i'^)W(iy\ij)PUW\jvM^\jv) (4) 

of the effective transfer matrix T, where the inequality (V', TV')/(V ', V) > 
(V,TV)/(V,V) is satisfied. 

Second, we decompose V into a matrix product 

v'ii'tW) - E A W\ a ) ^ AT Wv') (5) 

a 

via the singular- value decomposition (or the matrix diagonalization), where 
A is a nm by m' orthogonal matrix. We decide the new dimension m' for 
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the block-spin variable a (< m') so that the new singular values cu a in Eq.5 
are greater than a certain threshold. 

The third step is the renormalization of the effective transfer matrix 

£ AT ^W)W(i/j'\ij)P(j' V '\jv)A(jv\0 - P(*'m), (6) 
jj'vv' 

and that of the orthogonal matrix R 

Y, AT ^OR{K\p)A{ip\v) - Rim, (7) 

Kp 

where the graphical representation of these equations are shown in Fig. 2. 
Equation 7 is actually the renormalization of the product wave function, 
because the new R gives a new product wave function. After the renormal- 
ization by Eq.6 and Eq.7, P is a nm'-dimensional matrix, and R is a nm' by 
m! orthogonal matrix. At this point, we return to the first step by setting 
m = m' . 

We repeat the iteration shown above until R becomes invariant under 
the renormalization in Eq.7; this means A = R. At the fixed point of the 
DMRG method, the orthogonal matrix R satisfies the relation R(ia\f3) = 
u a R T '(a\i/3), where the matrix R(ia\f3) Up is proportional to the square of 
the corner transfer matrix. 9 - ) After R is converged to its fixed point value, 
we calculate physical quantities by using the fixed-point value of the product 
wave function V. 2,3 ^ 1 The matrix R gives correlation functions; Ostlund and 
Rommer 7 ) have shown that the largest eigenvalue of the matrix t{a^\^rj) = 
R(ia\j)R(i/3\r]) is equal to unity, and the second one gives the correlation 
length. 
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The convergence of R is rather slow when the correlation length is very 
long. In such a case, we modify Eq.4 to accelerate the convergence. In 
addition to V = TV in Eq.4, let us create another vector V" = TV. 
It is apparent that a certain linear-combination between V and V" gives 
a better Ritz vector. Such an improvement is called 'modified Lanczos 
method.' 12 ' 13 ) What we have to do is to substitute the best linear combi- 
nation V" = aV + bV" into the l.h.s. of Eq.5, where a and b are adjusted 
so that (V",TV")/(V", V") takes its maximum value. 

The PWFRG method and the DMRG method gives the same result at 
their common fixed point. The main difference between them is that the 
PWFRG method renormalizes the product wave function, while the DMRG 
method does not. With the use of the wave function renormalization (Eq.7), 
the PWFRG method avoid the numerical diagonalization of T. As a result, 
the PWFRG method runs much faster than the DMRG method. 

As a test case, we apply the PWFRG method to the square lattice Ising 
model. Figure 3 shows the nearest-neighbor spin correlation function of the 
Ising model. Numerical error in the data when m = 40 is less than 10~ 7 
outside the region 2.2 <T/J< 2.3, where J is the coupling constant. The 
result agrees with that obtained by the DMRG 6 ) and the CTMRG 9 ) methods. 

Finally, we discuss the way how to apply the PWFRG method to ID 
quantum lattice models. The transfer matrix of a 2D classical model is a 
product of Boltzmann factors, while the Hamiltonian of a ID quantum lat- 
tice model is a sum of local operators. Thus we can apply the PWFRG 
method to ID quantum models by replacing the transfer matrix renormal- 
ization in Eq.6 by the renormalization of the effective Hamiltonian. Such a 
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renormalization algorithm for the effective Hamiltonian is already given by 
the DMRG method 2 ' 3) . 

The authors would like to express their sincere thanks to Y. Akutsu and 
M. Kikuchi for valuable discussions. T. N. thank to S. R. White for helpful 
comments and discussions about the DMRG method. The present work is 
partially supported by a Grant-in-Aid from Ministry of Education, Science 
and Culture of Japan. Most of the trial numerical calculations were done by 
NEC SX-3/14R in computer center of Osaka university. 
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Figure Captions 

Fig. 1. Graphical representation of the effective transfer matrix T in Eq.2 and 

the product wave function V in Eq.3. 
Fig. 2. Graphical representation of the renormalization for P and R in Eq.6 and 

7, respectively. 

Fig. 3. Nearest neighbor spin correlation function of the Ising model. 
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